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ABSTRACT 


The statistical model assumes a statistical distribu¬ 
tion of nuclear density levels in the decay of the compound 
nucleus. The densities of nuclear levels can be described 
by a statistical formxila in a region containing many energy 
levels. This formula can be used to predict the cross- 
sections for photonuclear reactions within the energy 
range 0 to 30 MEV. Provided assumptions are made 
neglecting the coulomb effect, spin, parity, angular mom¬ 
entum, and tertiary and higher order reactions, these 
predictions can be made by a digital computer and the 
resvdts plotted on a graph with experimental data for 
comparison. The digital computer program is used to find 
the statistical level density formula that most closely 
approximates the cross-sections obtained experimentally. 
The formulas used thus far do not approximate the exper¬ 
imental data closely enough. Some formulas for the 
density functions give a shape for the cross-section as 
a fiinction of energy that corresponds well with the 
experimental curve, but the ordinates are too small. The 
formiila P(E) = C E EXP (A E^) is the best one found 


so far. 
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Introduction 


A description of a nuclear reaction in terms of the 
decay of long-lived equilibrium systems which neglect phase 
relations can be called a "statistical model". It is 
recognized that statistical theory qualitatively predicts 
the major portion of nuclear reaction cross sections and 
in this thesis a comparison is made between the statistical 
model and experiment. The reason for this is primarily 
that the decay of an equilibrium system is governed by 
the available phase space. ^ The use of nuclear models is 
an attempt to explain in tangible form the observed 
phenomena characteristic of nuclear reactions. 

To explain large resonances in neutron cross-sections, 
Niels Bohr in 193^ introduced the concept of a compound 
nucleus which is a many-body system of strongly interact¬ 
ing particles formed by the combination of an incident 
particle ^ with a target nucleus: 

X + a-*- C* 

In such a system the incident particle a has a short mean 
free path and shares its energy with the other particles 
of the system C* so that it cannot be re-emitted until, 

^Ericson, Torleif, The Statistical Model and Nuclear 
Level Densities. Advances in Physics, v. 9> no. 36 , 

Oct. i 960 : 430 . 
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as a result of fixrther exchanges, sufficient energy can 
again be concentrated on this or on a similar particle. 
Should the incident particle be a slow neutron, this could 
be a length of time sufficient to permit the emission of 
radiation. The capture process is complete, and the 
compound nucleus C is formed in its ground state either 
by emission of a single photon 

C* —► C + V 

or by the emission of a sequence of qxxanta in cascade. 
This explains the predominance of the capture reaction 

and the suppression of elastic scattering.^ 

The progress of many types of nuclear reactions 

can be described in terms of the compoiind nucleus as a 
two-stage process: 

X + ►Y + b 

in contrast to a single-stage process: 

X + a —►Y + b 3 

Niels Bohr observed the usefulness of dividing a 
nuclear reaction into two states: 

(1) the formation of the compound system 

(2) the disintegration of the compound system 
into the reaction products. 

^Blatt, John M., and Weiskopf, Victor F., 
Theoretical Nuclear Physics. John Wiley and Sons, New 
York. 1952: 528-529. 

^Blatt and Weiskopf: 528-529. 
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These two states can be treated as independent processes 
in the sense that the mode of disintegration of the com¬ 


pound system depends only on its energy, angxilar momentum, 
and parity, but not on the specific way in which it was 
formed. This assumption has proved to be valid or approx¬ 
imately valid for many cases.^ 

2. Theory of Nuclear Reactions 

a. Channels 

Consider the following nuclear reactions; 


a + X 


X + a (elastic scattering) 
X* + a 
Y + b 

,2 + c' etc. 


All of these reactions except that of elastic scattering 
can be subdivided according to the quantum state of the 
residual nucleus and the emerging particle. One can 
denote the states of these nuclei by o( ^ j ,-—, 

and the states of the incident or emerging particles by 

c'n" j j ,-,. Should particles a,b,c,-, be 

elementary, states ci("j 1^*' j -» refer to their 

spin orientation, and the following reactions result; 

X^i + a^ii 

Jr' + 




^Blatt and Weiskopf; 341* 


p 

Ji 
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(^* and o(^" are states of the target nucleus and the 
incident particle; thus,^* and or y'and can denote 

any quantum state of X and a or Y and b respectively 
which can occur in this reaction. The conservation laws 
of energy, angular momentum, and parity restrict the 
possible combinations of ^ ' and^",y' and y" , etc. 

Any such possible combination of residual nucleus and 
emerging particle, each of which is in a definite quantum 
state, is called a reaction channel. These channels are 

denoted by single Greek letters o(»^» Y~,-» which 

consist of both indices and or^’ and^". Channel 

X^i + a ^" is called the entrance channel or initiating 
channel of the above reaction.^ 
b. Reaction R'obability 

According to the Bohr assumption, one can write 
the cross-section of a nuclear reaction X(a,b)Y in the 
form: 

6'(a,b) =6'c(a) ^^(b) 

where 6^(a) is the cross-section for the formation of a 
compound system by particle a incident upon the target 
nucleus X. 6^(b) is the probability that the compoxind 

^Blatt and Weiskopf; 313 • 


k 





system C, once formed, decays by emission of a particle 
b, leaving a residual nucleus Y. (5^(b) is a pure number; 
the compound system C miist eventually decay in some 
special way; (?-^(b) is the probability for this special mode 
of decay.^ 

It is useful to describe the reaction X(a,b)Y in more 
detail by considering the cross-section corresponding 

to a specific entrance channel d\ and a specific exit 
channel ^ » Therefore, one specifies the quantum states 
of all reaction pairs before and after the reaction. Thus: 

e(d,/9) (^) where is 

the cross-section for the formation of C through channel 
, and is the probability of C decaying throtigh 

channel'^ 

The probability (^(/^) that the compound nucleus C 
will disintegrate throtigh a specific channel ^ can be written 
in the form 



-v 




k 


444 


(v-) 


where k = X and the sum is extended over all channels V 
into which the compound nucleus can decay. Thus the 


^Blatt and Weiskopf; 342. 
"^Blatt and Weiskopf: 342. 
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cross-section for the reaction can be calculated if 

the cross-sections ^(V)for the formation of the compoimd 
system by all possible channels are known and if the Bohr 
assumption is valid. The probability is energy dependent, 

and its value for a compound nucleus excited to an energy 
Ej, can be obtained from the above expression provided 
the cross-sections €^(y') and the wave numbers k<yr are 
known for those channel energies Sy "Wfhich produce the 

O 

compound nucleus with an energy E^. 

c. Decay of the Compound Nucleus 

The kinetic energy in channel ^ is given by 
C — C + /O where is the "Q" value from channel 
to channel p . Because most of this energy is the 
kinetic energy of the emitted particle b and only a small 
portion is the recoil energy of the residual nucleus Y, 
one can consider to be the kinetic energy of b, and; 

- i'cn - 

where £|jy Is given A - £bY Qab ^ 

is the kinetic energy in the entrance channel ( both a 
and X are in their ground states). This represents the 
maximum value of C which can occur if the residual 

T 

®Blatt and Weiskopf: 344-345 • 
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nucleus Y remains in the ground state. is the excita¬ 

tion energy of Y after the reaction and is the difference 
between the energy of Y after the reaction and the 
energy of the ground state of Y. So long as E^ is 
less than c^jy’ every level E^ of the residual nucleus 

there corresponds an energy(C of the emerging particle.*^ 

r 

Should the energy be high enough, many levels of the 
residual nucleus can be excited. The energy distribution 
of the particle b becomes continuous when the energies 
E^ are closer together than the definition of the energy 
in the incident beam or closer than the energy resolution 
of the experimental arrangement for detection of the 
reaction products. The shape of the distribution function 
®b(£) d£ for the particles b emitted between £ and £+d£ 
is given by: 


Gb(^:)d£ = IGc(/S) 

€ ^ €p< £ t <^.£ 

where the sum is extended over the channelsy^ whose 
energy lies within the energy interval d£ . The number 
of terms in this sum is given by the number of levels of 
the residual nucleus Y with an excitation energy E^ 
between E and E - d£ , where E “ ^b Y " ^ • One can 


9 


Blatt and Weiskopf: 


365 - 366 . 
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call this number 'V\^(E)d^ and Wy(E) the "level density". 
In order to determine the relative intensity distribution 


of the emerging particles I^(£')d£' » one need consider 


only the numerator in the expression for the probability 





^6^ ) = is a function of the channel 

The factor replaces the factor 


Experimentally, it is impossible to determine the 
channels by ■which a particle b left the nucleus or even to 
determine its energy. One is therefore interested primar¬ 
ily in the cross-section of an (a,b) reaction regardless of 
the specific channels. Should the energy of the incident 
particle be high enough to leave the residual nucleus in 
many different possible excited states, one can express 
the probability in terms of the level density W. Letting 



where the sum extends over all open 


channels^ leading to the emission of b, this sum can be 
expressed in terms of an integral: 



^*^Blatt and Weiskopf: 367 . 
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compound nucleus by a collision 'with the energy between 
b and the excited nucleus Y*, the latter having an excita¬ 
tion energy of<^Y These cross-sections depend 

only upon the channel energy^ « ^ function of the 

r 

maximum energy the particle b can acquire and 

leave the residual nucleus Y in its ground state. If the 
level densities Wy(£^) are known, the function F-j^ can 
be calciilated for all particle types. Thus the cross- 
section for an (a,b) reaction can be expressed: 

^(a,b) = ^(a) Fk 

i 

where the sum extends over all particle types c that may 
be emitted in this reaction. 

The average energy E of a system is a monotonically 
increasing fimction of the nuclear temperature 6. The 
fxmction E(0) has a vanishing derivative /dE \ =o for 6=o. 

If one expands E(e) near 0=o in powers of 0, it must 
begin at least with a quadratic term. This leads to the 
following assumptions: 

(1) E = E(0) permits a power series expansion 
arovind 0 = o. 

(2) the leading term of this expansion is propor- 
tional to 0 . 
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Blatt and Weiskopf: 369-370. 
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(3) the higher order terms of 0 can be neglected* 


One can now write, with a = constant, E = a 6*^ from 


which: 



0(E) 


= 2(aB) + constant)and for 


the level density: 


W(E) = C EXP(27S[^ ) 


This is a statistical formixla which applies to a region 
containing many energy levels. It is a very rough 
approximate orientation of the general features of the 
densities of nuclear levels. These level densities can be 
measured experimentally by comparing the energy distribu¬ 
tion of the emerging particles in a nuclear reaction with 
the results obtained from the theoretical formula. 

d. Secondary Nuclear Reactions 

In most nuclear reactions of the type a + X =Y + b, 
the residual nucleus Y is left in an excited state. The 
excess energy is emitted by Y many times in the form of 
one or more gamma-rays until such time that the nucleus 
reaches the ground state. In some instances the excita¬ 
tion of Y may be high enough for Y to emit another 
particle, perhaps a neutron or a proton. One can say 


^^Blatt and Weiskopf: 372. 
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then that the first nuclear reaction is followed by a 
••secondary nuclear reaction^': 

Y = S + c 

in which a particle c is emitted. The observed reaction 
is a bombardment of X by a, followed by an emission of 
b and c leaving a residual nucleus 2. One can represent 
this process as an X(ajbc)2 reaction. The most common 
reactions of this type are the (a}2n) reactions. 

Observed X(a}bc)3 reactions may be caused by two 


different mechanisms: 



X + a = Y + b 

Y = 2 + c 

(A) 

X + a = Y^ + c 

Y^ = 2 + b 

(B) 


If the time between the emission of the first and the 
second particle should become very short (of the order 
of the transit time of nucleons over nuclear dimensions), 
one can no longer refer to a defined excited state Y or 
Y^ between the two emissions. The reaction should then 
be described: 

X + a- 2 + b + c (C) 

This third process constitutes a simultaneous breakup of 
the compound system C into a nucleus 2 and two particles 
b and c. Process (C) does not correspond to any of the 

^^Blatt and Weiskopf; 374 • 
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"channels" which have been used so far. Blatt and 
Weiskopf indicate that an assumption should be made that 
the reaction X(ajbc)S does not go to process (C). There 
is some doubt as to whether or not this reaction should 
be included; however, in the computer calculations which 
follow, this assumption is not made and the process (C) 
is included. 

The cross-section for process (A) is denoted by 
^ (a;bc) and the cross-section for process (B) by 
^ (a;c,b). The observed cross-section of the (a;bc) 
reaction will be called §'(a;bc). It is the sum of the 
cross-sections for processes (A) and (B): 

^(a;bc) = 6‘(a;b,c) +^(a;c,b) 

Usually one of these two contributions is much more 
important than the other, but there may be cases in 
which 6'(a;b, c) and 6'(a;c,b) are of the same order of 
magnitude. The distinction between processes (A) and 
(B) should not be made if b and c are the same kind of 
particle (e.g., (a;nn) reactions). 

Considering the definition 

6'(a,b) = ^^(a)GQ(b) 

^^Blatt and Weiskopf: 374 - 375 • 

^^Blatt and Weiskopf: 375 * 
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the cross-section (a,b) also includes those reactions in 
which the residual nucleus undergoes a secondary reaction 
at a later time. One can therefore distinguish the cross- 
section 6*(a,b) which constitutes only those reactions 
where the residual nucleus does not emit additional 
particles, 6^(a,b) is what is vinderstood conventionally 
by the cross-section of an X(a,b)Y reaction. The follow¬ 
ing relation holds: 


6"(a,b) =^*(a,b) +£6-(a;b,c) 

c 

One can express the cross-section for the (a;b,c) 
reaction (process (A)) as follows: 

6'(a;b,c) =.^ (o() 

(X is the entrance channel, and the prime on the sum 
indicates that it is extended only over those channels ^ 
in which the residual nucleus Y is excited enoiigh to emit 
a particle c. is the relative probability (a pure 

number) that the compound nucleus C will decay into 
channel^, while Gy^'(c) is the relative probability that 
the nucleus Y in state^' (corresponding to channel p) 


will emit a particle c, 
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The simplifying assumption can be made that the 
residual nucleus Y always emits a particle if its excitation 


^^Blatt and Weiskopf: 376, 
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energy is high enovigh. Let be the smallest 

separation energy of a particle from Y. The assumption 
is then expressed by: 

^ provided that £^y “ ®min^^^* 

This is the eqxxivalent of neglecting the emission of gamma- 

rays in competition with the emission of particles, an 
assumption which is approximately correct in most cases 
except in the vicinity of the threshold. One cau thus 
write: 

^*(a,b) 

where the double prime on the sum indicates that it is 
extended only over those channels in which the residual 
nucleus Y is excited below Y). Let ^sec 

maximum energy available for emission of secondary 
particles from the first residxxal nucleus Y for a given 
excitation energy in the original compovind nucleus. This 
energy is also the energy of the emerging particles b in 
the primary X(a,b)Y reaction at which Y can just begin 
emitting particles (See Figure 1). ^ is given by: 

S 0C 

^sec ~ ^bY “ ®min^^^ 

One then obtains the following expressions for the 
"observed” cross-section g'^Cajb) of the X(a,b)Y reaction: 
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Energy Relationships in the Reaction 
X -h a z Y* + b , 

!♦ = Z -h c 

FIGURE 
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^*(a,b) = 6'(a,b) 


6'*(a,b) 


4 


(for £ < o ) 

^sec 


:e(a 


cf Ib 

,b) ^sec (for£‘>o) 



Where 1^(6^) d^" is the relative probability of emitting 
particles b with energies X between £ and (£ + d£ . 

( = constant£^(^) ^ 

The cross-section for a secondary nuclear reaction 
with emission of particles c following a primary (a,n) 
reaction can be written with the same approximation as: 


£'(a;n,c) ^ C(a,n)_o 


exp( - )Gy/ 9 ' (c)d(f 


:bY 


exp( - £ )d£ 

e 


= fbY - (Y) 

In this formiila is the excess energy over the threg^. 

20 

hold of the (a;n,n) reaction. S^(Y) = Sj^(Y). 

Shoxild the energy of the incident particle be 
svtfficiently large, the residual nucleus after emission of 
the second particle (nucleus 2) may be excited eno\igh to 


^"^Blatt and Weiskopf: 376. 
^®Blatt and Weiskopf: 377* 
^^Blatt and Weiskopf: 377. 

^^Blatt and Weiskopf: 378-379. 
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emit a third particle. One would then observe tertiary 

21 

reactions or reactions of even higher order. 

3. The Statistical Model 

The model portrayed herein is a statistical representa 
tion of the nucleus and has been adapted for calculations 
in a digital computer. This model, with due regard to 
the appropriate assumptions, is general and can be used 
to make statistical predictions regarding the cross-section 
for the emission of particles from any initial nucleus. In 
its current status, the computer program is somewhat 
limited in that the following nuclear properties have been 
neglected: 

a. The coulomb barrier 

b. Spin 

c. Parity 

d. Angular momentum 

e. Tertiary and higher order reactions 
The density function P(E) is the same for the 

different product nuclei. 

Tertiary reactions were considered, but have not 
been included in the main program. There is included in 
Appendix II a computer subroutine which can be used to 

^^Blatt and Weiskopf: 378 
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compute tertiary reactions, 

I 

The coulomb force can affect the neutron processes 
indirectly. Should a proton process be made smaller, then 
the neutron process will be made larger. In the calcula¬ 
tions studied herein, the coulomb effects are neglected. 

In setting up the model and fitting it to the computer 
program, Cu^^ was selected as the experimental nucleus 
with which to compare the theoretical results. Cu^^ 
was selected for the following reasons: 

a. Cu^5 is a light nucleus and is less complicated 
than others, 

b. Cu^^ has a low coulomb barrier, 

c. Experimental data is available for making 
desired comparisons with theoretical calculations. 

The experimental data was taken from graphs of 
the published results of experiments with photo neutron 
cross-sections conducted at the Ernest O, Lawrence 
Radiation Laboratory, Livermore, California in August 
1963. (See Figures 2,3 and 4) In the Livermore experi¬ 
ments the single neutron cross-sections actually consist 
of the sum of the reactions: 

^np ^n2P + — — — 
and the double neutron reactions: 

^2n ^2np ^2n2p - 
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The general technique here is to compute the photo 
neutron cross-section for the X(y,2n )Y process and 
compare it with the experimental data for the same 
process. This is an application of the relationship: 



^(a;n,c)^6'(a,n) 0 




^(a,n) is taken from the experimental data for Cu^^ 


and is used throughout in all computer calciilations. For 
these processes, 6‘(a,n) is the and 6(a}n,c) is 

^(■v^;2n). 

The range of gamma-ray energy for which calcula¬ 
tions are made extends from 0 to 30 MEV. Within this 
energy range, six reactions are possible. These reactions 
are portrayed in Figure 5* 

In order to make computer calculations less complicated, 
channel widths are set at 0.3 MEV, and all integrals are 
calculated in accordance with this incremental change in 
energy. Also, the threshold energies for each reaction 
have been adjusted so as to be divisible by 0.3 for ease 
in carrying out the necessary numerical integration. 

The probability of decay to any particular state is 
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statistical, and is computed for both single-step and two- 


step processes. In setting up the model, the return of 
the compound nucleus to its ground state by the emission 
of a gamma-ray is neglected. 

From Figure 5» it is noted that there are two 
direct processes by which two neutrons can be emitted: 

a. Cu » Gu^ 3 + 2n 

b. Cu^5_^Ni^^ + p + 2n 

For simplicity, the Cu^5 ( 7 ';p, 2 n) Ni^^ process will be 
omitted. These omitted processes contribute to the 
overall cross-section only at energies above their thres¬ 
holds. For this process, the threshold is taken to be 
23.7 MEV. 

Four secondary processes res\alt in the emission of 


two neutrons: 

a. Cu^^ 

-►Cu^^ 

+ 

n—►Cu^^ + 2n 

b. 

Cu^^ 


+ 

p—►Ni^^ + p + 2n 

C 0 

Cu^5 

^Ni^^ 

+ 

p + n-*-Ni^^ + P + 

d* 

Cu^^ 

-►Cu 

+ 

2n-KNi^^ + p + 2n 


Two tertiary processes result in the emission of 
two neutrons: 

a. Cu^-5_>Ni^^ + p_>.Ni^3 + p + p 

b. Cu^-^->Cu^^ + n^-Cu^^ +2n->Ni^^ + P + 


+ 2n 
2n 
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In each case, the tertiary processes and the 
(Tjp,2n)Ni^^ processes are omitted since they contribute 
to the overall cross-section only at energies above the 
threshold of 23.7 MEV. 

Single neutron emitting processes are as follows: 
Three direct 


a. 

Cu^-^ 

-►Cu^^ 


n 

b. 

Cu^^ 

-►Ni^3 

+ 

p + n 

c. 

Cu^5 

-►Co^^ 

+ 

2 p + n 


The ('Vi2p,n) process is omitted. 
Three secondary 


a. 

Cu^^ 


+ 

P — 

-►Ni^^ + 

P 

+ 

n 

b. 

Cu^^ 

-►Ni^^ 

+ 

P — 

-►Co^^ + 

2 p 


n 

c. 

Cu^5— 

-►Cu^^ 

+ 

n — 

-►Ni^^ + 

P 


n 


The (VjSpjn) process is omitted. All tertiary single 

neutron processes are omitted. 

The decay of the compound nucleus for each process 

may be expressed as an integral of the density function 

integrated over all available energy channels: 

^E-Ethl 

j P(E')dE» 

0 

This integral represents the decay of Cu^^ to Ni^^ + p. 
Other reactions are represented accordingly. 
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Probabilities for direct processes are calculated as 


follows; 


a. Two-neutron processes 


(1) Cu^ 


“Cu 


,63 


2n 

E-Eth4 

...3 ) P, 

E-^:thl E-Eth2 E-Eth3 E-EthU E-Eth5 E-Eth6 
'p(E'clE'+ /p(E')clE'+/p(E')dE'+/p(E')cffi'+|P(E')dE' + /p(E')dE' 
0 ^0 ^0 ^0 ^0 ^0 


Where Pj(E) is the probability of decay of the compo\ind 
nucleus C* to Cu^^ with the emission of two neutrons. 

The function R( E - Eth5 ) will be explained in subsequent 
discussion. 


(2) Cu 


65_ 


■Ni + p + 2n 

E-Eth5 



P3(E) = 


(1) Cu^^-+ n 

E-Eth2 

/ R( E - Eth5) P(E')dE' 

_ 

E-Ethl E-Eth2 E-Eth3 E-Ethl; E-Eth5 E-Eth6 

P(E OdE' Y P(EOdE' y^(EMdE'+/^' )cffi' 4 /p^EO dE' h/p(E' )dE' 
0 0 ' 0 ' 0 ' 0 ^ 0 
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(2) Cu^-^ 



+ p + n 
E-Eth3 


P^(E) = 



R( E - Eth5) P(E') dE 



E-Ethl E-Eth2 E-Eth3 


E-Eth3 


E-Ethii 



E-Eth6 


P(E' )dE'+/p(E' )dE«+/P(E' )dE'-/P(E' )dE‘+/p(E‘ )dE'+/P(E' )dE‘ 


(3) Cu^- 



+ 2p + n (omitted) 


The probabilities for secondary processes are 
calculated with a slightly different procedure. These 
two-step processes consist of a decay of the compound 
nucleus to a second excited nucleus. This second excited 
nucleus then decays by emitting a particle. Thus the 
probability for the complete process will be formed by the 
product of the probabilities of two independent reactions. 
For the general scheme of decay, consider a hypothetical 
example. (See Figure 6) 

The available energy of the initial reaction ranges 
between 0 and 30 MEV and is represented by E. Decay 
of the compound nucleus to state A can occur if E is 
greater than the threshold energy EtA. Likewise, if E 
is greater than EtB, the compoxmd nucleus may also 
decay to state B. Should this condition prevail, then 
both states A and B will compete for the reaction. 
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Probabilities for the two-step process must then make 
provisions for the consideration of this competition. Now, 
once the initial reaction has taken place and a second 
compovind nucleus has been formed, this second compound 
nucleus will tend to decay if the energy is high enough. 
Thus, if E - EtA is greater than EtB - EtA, this 
second compound nucleus will tend to decay. The decay 
of this second compound nucleus by the emission of a 
particle will cause a decrease in the probability of this 
second state existing when the energy exceeds the value 
of E - EtB. This "leakage” must be compensated for 
when the probabilities are calculated. 

In calculating the probabilities for two-step reactions, 
the reaction which occurs second physically must be 
computed first mathematically. To compensate for the 
"leakage" factor, the function R( E - Eth) is defined 
as follows: 


R( 

E - Eth) = 0 

for 

E - Eth > 0 

R( 

E - Eth) = 1 

f or 

E - Eth <0 


The particular value of Eth will depend upon the reaction. 
The function R( E - Eth) becomes a part of the density 
function for the second reaction and is integrated along 
with the density function. 
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The "leakage” must also be considered for the single 


neutron and single proton processes for energy values 
above the two—neutron threshold since the two-neutron 
process will "rob" the single neutron process. This 
"leakage" is compensated for in the appropriate eqviations 
by the same fxmction R( E - Eth). 

Consider the two-step reaction; 


+ a + b 


+ a 


Where X represents the initial compound nucleus, Y 
represents state A, and 3 represents state B. Once 
state B has been formed, this state will exist xintil the 
energy E2 exceeds the value EtC-EtB, at which energy 
there will tend to be a "leakage" of particles from state 
B to state C (provided of course that state C is one 
which is physically possible from state B). This "leakage" 
is compensated for by the fiinction R( E - Eth). Thus 
the integral for the second of the two processes would 
be: 


E - EtA - (EtB - EtA) 

R( El+EtA - EtB) P(E')dE' 



0 


But, E - EtA -(EtB - EtA) = E - EtB, and El + EtA 
= E. Therefore: 


E - EtB 
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The integral for the complete two-step process, combined 
with the appropriate denominator, would appear thusly: 


E - EtA 


E - EtB 






>^^(E1') dEl» dE' 


0 


0 


D1 


D2 


D 1 represents the denominator for the decay of the 
compound nucleus to state A, and D2 represents the 
denominator for the decay of state A to state B. 

Using the above procediare, the probabilities for the 
two-step processes can be written. Consider first the 
secondary reactions: 


a 


ISii^if* + p-►Ni^^ + 


p + n 


R( E - Eth5) P(El') dEl' 





+ n + p 


R( E - Eth5) P(E2’) dE2' 



E - Eth3 E - Eth5 E - Eth6 




0 


P(E2')dE2 (E2')dE2' (E2’)dE2’ 



0 



0 
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Now, the probabilities for the two-step single-neutron 
processes can be written: 


a. Cu 




-Ni^^* + p —►Ni^^ + p + n 

E - Ethl E - Eth3 



P2pn(El) dEl dE' 

Po__(E)= ^0 ^0 

i:-Eth2 fc-Eth3 &-EthU E^th^ E-Eth^ 

P(E' )dE'Y P(E')clE»+/p(E')dE'+/p(E')ciE‘+/p(E' )dE'+/p(E' )dE' 

0 -^0 •'0 "0 *^0 '0 

b. Cu^5*-►cu^^* + n—►Ni^^ + n + p 

E - Eth2 E - Eth3 
P( E') 

0 ' 0 


P^ (E2) dE2 dE* 
2np 


^2np(^) = 


E-Ethl E-Eth2 E-Eth3 E-Ethii E-Eth5 E-Eth6 
P(E' )dE'+/p(E‘ )dE'+/p(E' )dE'H/p(E')dE'-/p(E' )dE'+fp{E* )dE' 

0 -'0 -'O *^0 *'0 ''0 

The probability for the two-step two-neutron 
process can be obtained similarly. Again consider first 
the secondary reaction: 

Cu^^* + n—►Cu^^ + 2n 

^2nn(E2)= R ( E - Eth5) P(E2')dE2' 


E-Eth3 E-Eth4 E-Eth5 

P( E2')dE2'+jP( E2')dE2'+/p(E2')dE2 


E-Eth6 

2^yP(E2')dE2' 
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The probability for the complete two-step process 


becomes: 

Cu^5*—^ Cu^^* + n—►Cu^^ + 2n 

E - Eth2 E - Ethlf 

/ P(E') / 

Po__(E)= "^0 >^0 

El-Ethl E-EthiJ E-Eth3 E-EthU E-Eth^ E-Eth6 

/p(E' )dE' + /T(E' )dE'+/p(E' )dE'+/P(E‘ )dE'+/p(E' )dE' + ^(E« )dE' 

*^0 ^0 ^0 ^0 ^0 ^0 


The total probability for each process is the sum of 
the probabilities for the individual reactions since these 
are mutvially exclusive events. The total probability for 
the two-neutron processes can be written: 

^nn<^l ' ^2nn(B) + ^i(E) + P2(E) 

The total probability for all single-neutron processes is: 

Pn(E) = P2p„(E) + P2np(E) + P 3 (E) + P^(E) 
Statistically, the ratio of the cross-sections for 
the two—neutron processes to the cross-sections for the 
single-neutron processes should equal the ratio of corres¬ 
ponding probabilities: 


and 


^ (V,2n) ^ P(y.2n) 

^ (Vfn) P(V,n) 


^(r,2n) =^(V,n) 


P(V,2n) 

P(V,n) 
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where P(y',2n) = ^(T»n) - Pjj^(E). This 

relationship is the basis on which the computer program 
is written. 

4 . Computational Parameters 

Blatt and Weiskopf present the function 

P(E) = C expCS^AE ) 

as a first guess in the determination of the density 
f\inction describing the statistical decay of the compound 
nucleus. In this function, E represents the energy of 
the incident particle (gamma—ray); C is a constant of 
proportionality; and A is a constant based on the nuclear 
temperature 0. 

From this first guess, several variations have been 
devised in an attempt to obtain a density function which 
will closely approximate the available experimental data 
for the (V,2n) process. In each of these density functions 
the following constants are used: 

A = Nuclear temperature constant; the units of 
A are (MEV)”^ 

C = Constant of proportionality 

^^Fultz, S. C., Bramblett, R.L., Caldwell, C. T., 
and Harvey, R.R., Photo-Neutron Cross-sections for 
Natural Cu, Cu^3, and Cu°5. Physical Review, v. 133» 
no. 5B,9 March I 964 : BII 49 - BII 54 . 

^^Fultz, Bramblett, Caldwell and Harvey:B1149-B1154* 

^^Fultz, Bramblett, Caldwell and Harvey:BII 24 . 9 -BII 54 • 
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n = RN = Exponent of the (AE) term 

Altogether, fourteen density ftinctions have been 
formvdated for this program (See Appendix III)j however, 
only six have been investigated in any detail. These six 
are? 

a. P(E) = C(3 - E) EXP (A( E^ - ) ; 3 = E-Eth 

b. P(E) = C EXF( A(E^)) 

c. F(E) = CS EXP( A( E’^)) 

d. P(E) = CE EXP( A( E’^)) 

e. P(E) = CE/EXP(A ( E^)) 

f. P(E) = C(E)3 EXP(A (E^)) 

Each of these density functions has been tested 
with various values of A. In all cases, the parameter 
n = RN was set equal to 1.0. The values of A ranged 
between 0.0001 and 10.0 (MEV)”^. 

For each of the values of A, the cross-sections for 
the theoretical (V,2n) process were plotted versus the 
energy of the incident gamma—ray. On the same graph in 
each case were plotted the experimental cross-sections for 
the corresponding values of gamma-ray energies. The 
procedure used is to compare the theoretical cross-sections 
with the experimental values. The objective here is to 
find a density function P(E) that will give a close 




) 






approximation to the experimentally determined cross- 
sections. 

5. Results 

a. P(E) = C(S - E) EXP( A( E^-S^); 3=E -Eth. 
See Figures 7 through 14. The idea for this function 
was obtained from Ericson’s article on The Statistical 
Model and Nuclear Level Densities. At values of A 
ranging from 0.0001 to 1.0 (MEV)”^, it shows a shape 
that in general conforms to the shape of the experimental 
curve; however, for larger values of A the cxxrve develops 
a spike near the lower energy values. Thxis, this function 
is not considered to be appropriate. 

b. P(E) = C EXP(A(En)). See Figxxres 15 
through 23 . This fionction is a modification of the one 
suggested by Blatt and Weiskopf as a first guess. This 
fxinction was tested for values of A ranging from O.OOO 5 
to 10.0 (MEV)”^. This function reached its maximum 
value at an A value of about 0.1 (MEV)”^ and held this 
value for all decreasing values of A. As the value of A 
increased from 0.1 (MEV)”^, the fionction decreased in 
value. The values became so small at an A value of about 
2.0 (MEV)”^ that the curve did not even show up. In 

all instances, a small spike appeared on the lower energy 
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end of the curvej however, it was not as prono\uiced as 
the spike obtained with the ftmction in a. above. This 
function is not considered to be descriptive of the 
experimental cross-section. 

c. E{E) = CS EXF(A Ell), see Figures 24 
throiigh 33. For this function, A ranged in value from 
0.0001 to 10.0 (MEV)“^. Throvighout the entire range of 

A values, the shape of the curve was generally in conform¬ 
ity with the experimental c^lrve except that the ordinates 
were at all times much too small. This function reached 
its maximum value at an A value of about 0.5 (MEV)"^. 

For increasing values, the curve became ^lnobservable at an 
A of about 1.5 (MEV)”^» This fiinction gives a good 
shape, but the values of the cross-section are too small. 

d. E{E) = CE EXP( A( E") ). See Fig^lres 34 
through 49• Values of A ranged from 0.0001 to 10.0 
(MEV)”^. The shape of the curve was good throughout, 
but its maximum value never exceeded about one—fourth 
of the experimentally determined values. It reached its 
maximum value at an A of about 0.01 (MEV)“^J it became 
too small to observe at an A of about 1.5 (MEV)”^» 

Tj'his fTonction does not have a high enough maximum value. 
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e. P(E) = CE/EXP( A(E^)). See Figures 50 
through 56 . The values of A ranged from 0.0001 to 
10.0 (MEV)“^ for this function. Throughout this range 
of values, the curve was visible and had significant 
magnitude. However, there developed a spike on the low 
energy end which became quite pronounced as the values 
of A increased. This is not considered to be an approp¬ 
riate density function. 

f. P(E) = C(E)3 EXP( A( E^)). See Figures 57 
through 64 . This function was devised in an attempt to 
increase the magnitude of the ordinates for the cross- 
section curve. From the figvires, it can be seen that 
instead of increasing the magnitudes of the ordinates, it 
actually decreased them. For this function, values of A 
ranged from 0.0001 to 3.0 (MEV)”^. Throughout most of 
this interval, the curve held a good shape. The magnitudes 
became so small that as the values of A increased through 
1.0 (MEV)“^ the curve was no longer visible. Although 

the energy coefficient (E)3 does not adequately describe 
the process, it does suggest the possibility that perhaps 
some other power of (E) may give the curve a better 
shape. 
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Conclusions 


None of the density fvinctions tested thus far give 
an adequate approximation to the experimentally obtained 
data. Two of the tested functions produce a good shape, 
but in all cases the maximum values do not approach 
the experimental values. The lack of conformity so far 
may be due to several causes; 

a. The proper density function has not been found. 

b. Co\ilomb and other effects have been ignored. 

c. Tertiary reactions have been omitted. 

The results indicate that there may be some power 
function of (E) which acts as a variable coefficient of 
the exponential term. The density function with the (E)^ 
coefficient is a first attempt to find such a term. 

The results obtained, although not in close agreement 
with experimental data, are encouragng and indicate that 
this model with appropriate modifications and development 
can bs used for making useful statistical approximations 
for mclear reactions. 
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APPENDIX I 


CROSS-SECTION PROGRAM 
This program was written for the purpose of 
calculating cross-sections in the decay of the compotmd 
nucleus. Although written for the X(>^,2n)Y reaction 
with Cu^^ as the target nucleus, this program can be 
used for any nucleus in computing an approximation to the 
cross-sections for the emission of particles. Appropriate 
assumptions must be made with regard to the coulomb 
barrier, spin, angular momentum, etc. 

In general, the program consists of a series of 
calcidations in which three subroutines perform most of 
the mathematical operations. These subroutines calculate 
values of single and double integrals. A subroutine for 
the computation of triple integrals is included, but was 
not vised in the main program. 

Subroutine PROD calculates the values of single 
integrals. The call statement is: 

CALL PROD( EL, PEN, PNO, NWB) 

The variables in the call statement have the following 
meaning: 

EL is the threshold for the reaction being calculated. 

PIN is a constant used in calculating double and 
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triple integrals. For the single integral, PIN equals 1.0. 

NWB is the index corresponding to the energy range 
over •which the calculation is made. It is determined by 
dividing the energy in MEV by 0.3* (©•§•» NWB = E/.3). 

It is appropriate to note here that all threshold 
energies have been adjusted so as to be divisible by 0.3 
in order to simplify the computer operations. 

Subroutine GRASP calculates the values of double 
integrals. The call statement iss 

CALL GRASP ( ELI, EL2, EL3, EL4, EL5, EL6, 
PN20,POCN,PNLN,EFU,NPT,NP1,NP2,NP3,NP4,NP5, 
NP6) . 

In order to observe more clearly the meaning of each 
variable, consider a two step reaction: 

X(r,n)Y* ; Y*—►S + n 
For the total reaction of X +Y —►S + 2n over the 
energy range of 0 to 30 MEV, six reactions are possible. 
Let the probability for the first reaction be: 


E-Ethl 




E-EfEl 


E-Etli2 E-hth^ h-Ethh £-Eth5 E-Eth6 


0 
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No-w PNLN = PNLl corresponding to the Ethl for the 
threshold one reaction in the numerator. ELI through EL6 
are the threshold energies for each reaction that can take 
place over the prescribed energy range. EPL is the 
threshold energy for the initial reaction; in this example 
EHL = ELI. 

Next, consider the second reaction: 

E - Eth3 



In the call statement, POCN = POC3 corresponding to 
the Eth3 for the threshold three reaction in the numerator. 
EPU is the threshold for the second reaction and in this 
example EPU= EL3. The variables NPl through NP6 are 
indices for the computation of the integrals in the denom¬ 
inator. If a particular integral is not to appear in the 
denominator, then the corresponding index will be set 
equal to 2; if the integral is to appear in the denominator, 
the index will be 1. For example, in the probability for 
the second reaction the integral for Ethl does not appear 
in the denominator; therefore, NPl will equal 2 in the call 
statement. For the above probability, the sequence NPl 
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through NP6 'will appear in the call statement as 2,2,1,1,2,1. 

The value for NPT is the "leakage" factor. In this 
example NPT = Eth5/0.3« 

PN20 is the output of this subroutine and it 
represents the probability of the formation of cross- 
sections by a two-step process. PN20 may have any 
name in the call statement. 

The triple integral subroutine follows the same 
general procedxire. 
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APPENDIX II 
PROGRAM FLOW CHART 


START 


READ SEA 

e 

FOR ALL n 
PROCESSES 


READ SEC 
(5* 

FOR ALL 2n 
PROCESSES 


E(l) =0.0 
J = 1 


E(J-i-l) =E(J) + 0.3 


ML 
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APPENDIX IIB 


PROGRAM LISTING 


PROGRAM PRODUCE 


DIMENSION SEA ( 101) ♦SECdOl) , Y ( 102 ) , 1 T I T LE ( 12 ) , 

1 E(102) ♦PNU(101)tPNAI101)♦PNSI101) ,PNJ(101)» 

2 PNCKlOl) ,PNC2(101) ,PNC3 ( 101) ♦ PNC4 (101 ) , 

3 PNC5( 101) ,PNC6( 101 ) ♦PN7( 101 ) ♦PNG! 101.) jCRNI 101 ) ♦ 

4 PNR(lOl),PNK(101)»PNP(101),DENOM(101) 

COMMON E»P0C1.P0C2,POC3»POC4»POC5»POC6, 

1 PNL1,PNL2»PNL3»PNL4,PNL5»PNL6»RN,C,A» 

2 DENOM 
RN = 1.0 
C = 1.0 
A = 1.0 
READ 99 SEA 
READ 99 SEC 

99 FORMAT ( 20F4.0) 


Ed) 

= 0.0 


\ 


DO 5C 

1 J = 

1,101 



E(J + 

1 ) 

= E(J) + .3 



NWB = 

101 




PIN = 

‘ 1.0 




EL = 

7.2 




CALL 

PROD 

( EL,PIN,PNC1, 

NWB 

) 

EL = 

9.9 




CALL 

PROD 

( EL,PIN,PNC2, 

NWB 

) 

EL = 

17.1 




CALL 

PROD 

( EL,PIN,PNC3, 

NWB 

) 

EL = 

17.7 




CALL 

PROD 

( EL,PIN,PNC4, 

NWB 

) 

EL = 

23.7 




CALL 

PROD 

( EL,PIN,PNC5, 

NWB 

) 

EL = 

26.7 




CALL 

PROD 

( EL,PIN,PNC6, 

NWB 

) 

DO 6161 J 

= 1,101 



PNC22(J) = 

PNC2(J) 



PNC33(J) = 

PNC3(J) 



PNC44(J) = 

PNC4(J) 




6161 CONTINUE 

DO 990 J = 1»101 

DENOM(J) = PNCl(J) + PNC2(J) + PNC3(J) + 
1 PNC4(J) + PNCSIJ) + PNC6(J) 

IF( J - 57 >52,49,51 
49 I = J 

51 PNCX = PNC22(I) 

GO TO 57 

52 PNCX = PNC2(J) 
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57 PN6(J) = PNCX/DENOM(J) 

IF( J - 79 )59,53,54 

53 K = J 

54 PNCY = PNC33(K) 

GQ TO 67 

59 PNCY = PNC3(J) 

67 PN7<J) = PNCY/DENOM(J) 

IF( J - 79 )103,101.102 

101 M= J 

102 PNCZ = PNC44(.M) 

GO TO 77 • 

103 PNCZ = PNC4(J) 

77 PN4(J) = PNCZ/DENOM(J) 

990 CONTINUE 

CALL GRASP ( 7.2,9.9,17.1,17.7,23.7,26.7, 

1 PNP,POC3,PNL2,17.1,9.9,79, 

2 2,2,1,2,1,1 ) 

CALL GRASP ( 7.2,9.9,17.1,17.7,23.7,26.7, 

1 PNR,P0C3,PNL1,17.1,7.2,79, 

2 2 , 2 , 1 , 2 , 1,1 ) 

CALL GRASP ( 7.2,9.9,17.1,17.7,23.7,26.7, 

1 PNU,P0C4,PNL2,17.7,9.9,79, 

2 2 , 2 , 1 , 1 , 1,1 ) 

DO 1000 J = 1,101 

PNJ(J) = PN4(J) + PNU(J) ' 

PNK(J) = PN6(J) + PN7(J) + PNR(J) + PNP(J) 

CRN(J) = SEA(J) *(PNJ(J))/ PNK<J) 

1000 CONTINUE 

NUMPTS = .101 
MODCURV = 1 
Yd) = 100.0 
D041J=2,101 

41 Y(J) = CRN(J) 

ITYPE = 0 
LABEL = 4H 
ITITLE(l) = 8HL0WRY 
ITITLE(2) = 8HA = 1.0 
ITITLE(3) = 8HRN= 1. 

ITITLE(4) = 8HE=FNERGY 
ITITLEI5) = 8H1 UNIT 
ITITLE(6) = 8H= 1 MEV 
ITITLE(7) = 8HY=CR0SS 
ITITLE<8) =8HSECTI0N 
ITITLE(9) = 8H1 UNIT 
ITITLEI10) =8H= 10 MB 
ITITLEI11) = 8H F - 1 
IT ITLE(12) = 8H CU 65 

CALL DRAW ( NUMPTS,E,Y.MODCURV,I TYPE,LABEL,ITITLE, 

1 EXSCALE, YSCALE, IXUP, IRIGHT, MODEXAX, MODEYAX, IWIDE, 

2 IHIGH, IGRID, LAST ) 

MODCURV = 3 

DO 42 J =*2,101 

42 Y{J) = SEC(J) 

CALL DRAW ( NUMPTS,E,Y.MODCURV,I TYPE,LABEL,ITITLE, 
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1 EXSCALEt YSCALE* IXUP* IRIGHT. MODEXAX» MODEYAX* IWIDE* 

2 IHIGH» IGRID* LAST ) 

END 

SUBROUTINE GRASP ( EL1*EL2,EL3» 

1 EL4»EL5»EL6»PN20,P0CN»PNLN*EPU,EPLiNPt♦ 

2 NP1,NP2*NP3»NP4»NP5»NP6 ) 

DIMENSION E(102).PN20(101),DENOM(101) 

COMMON E*P0C1»POC2»POC3*POC4.POC5»POC6, 

1 PNL1»PNL2*PNL3»PNL4*PNL5*PNL6»RN,C»A» 

2 DENOM 

DO 889 J = 1.101 
IF(E(J) - EPL )10,3,3 

3 CONTINUE 

NWB = 10.*((E(J) - EPL)/3.) + 1.0 
EP = E(J) - EPU 
IF(EP)11.4.4 

4 CONTINUE 
PIN = 1.0 

1 EL = ELI - EPL 
IF(EL)22.12.12 

12 GO TO ( 121.22 ).NPl 

121 CALL PROBE II ( EL.PIN.POC1.NWB ) 

GO TO 2 
22 POCl = 0.0 

2 EL = EL2 - EPL 
IF(EL)28.13.13 

13 GO TO ( 131.28 ) .NP2 

131 CALL PROBE II ( EL.PIN.POC2.NWB ) 

GO TO 8 

28 POC2 = 0.0 

8 EL = EL3 - EPL 
IF(EL)25.14,14 

14 GO TO ( 141.25 ).NP3 

141 CALL PROBE II ( EL.PIN.POC3.NWB ) 

GO TO 5 

25 POC3 = 0.0 

5 EL = EL4 - EPL 
IF(EL)26.15.15 

15 GO TO ( 151.26 ).NP4 

151 CALL PROBE II ( EL.PIN .POC4.NWB ) 

GO TO 6 

26 POC4 = 0.0 

6 EL = EL5 - EPL 
IF(EL)27.16.16 

16 GO TO ( 161.27 ).NP5 

161 CALL PROBE II ( EL.PIN.POC5.NWB ) 

GO TO 7 

27 POC5 = 0.0 

7 EL = EL6 - EPL 
IF(EL) 29.17.17 

17 GO TO ( 171.29 ) .NP6 

171 CALL PROBE II ( EL.PIN .P0C6.NWB ) 

GO TO 30 

29 P0C6 = 0.0 

12j.O 


( 

I 

I 


( 




I 



II 

\ 


I 


( 


I 








30 CONTINUE 
GO TO 779 
11 PNA = 0.0 
GO TO 002 
779 CONTINUE 

IF( J - NPT ) 0OO»2O1,2O2 

201 POCP = POCN 

202 POCN = POCP 

000 PNA = POCN/ ( POCl + POC2 + POC3 + P0C4 + 
1 P0C5 + POC6 ) 

002 CONTINUE 
EL = ELI 


CALL 

PROBE 

I I 

{ 

EL,PNA,PNL1,J 

) 

EL = 

EL2 





CALL 

PROBE 

I I 

{ 

EL»PNA»PNL2»J 

) 

EL = 

EL3 





CALL 

PROBE 

I I 

( 

EL»PNA,PNL3»J 

) 

EL = 

EL4 





CALL 

PROBE 

I I 

( 

EL»PNA»PNL4»J 

) 

EL = 

ELS 





CALL 

PROBE 

I I 

{ 

EL»PNA»PNL5»J 

) 

EL = 

EL6 





CALL 

PROBE 

I I 

( 

EL.PNA,PNL6»J 

) 


PN20(J) = PNLN/DENOM(J) 
GO TO 001 
10 PN20(J) = 0.0 
001 CONTINUE 


POCl 

= 

0.0 

POC2 

= 

0.0 

POC3 

= 

0.0 

P0C4 

= 

0.0 

POC5 

a 

0.0 

POC6 

a 

0.0 

PNLl 

= 

0.0 

PNL2 

= 

0.0 

PNL3 

= 

0.0 

PNL4 

a 

0.0 

PNL5 

= 

0.0 

PNL6 

= 

0.0 

009 CONTINUE 


RETURN 

END 

SUBROUTINE PROD ( ELfPINtPNO.NWB ) 

DIMENSION E(102).P(101).PNOI101),DENOM(101) 
COMMON E» POC1»POC2»POC3»POC4,POC5»POC6» 

1 PNL1»PNL2»PNL3»PNL4»PNL5»PNL6»RN.C,A» 

2 DENOM 

DO '2 K =1»NWB 
Z = E(K) - EL 
IF(Z)10,3»3 
3 H = 0.3 

W = 10.0*{Z/3.) 

6 NW * W + 1.0 
Y * 0.0 
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DO 206 L =1*NW 

P(L) = C*(Z - Y )*PIN*EXPF(A*(Y**RN - Z**RN )) 

206 Y = Y + H 
SUM I = 0.0 
NY = W 

DO 207 L.=2» NY 

207 SUM I = SUM I + P(L) 

PNO(K) = H/2.*(P(1) + 2.*SUM I + P(NW)) 

GO TO 2 

10 PNO(K) » 0.0 

2 CONTINUE 
RETURN 
END 

SUBROUTINE PROBE II ( EL.PIN,PNO»NWB ) 
DIMENSION E(102).PI 101),OENOM(101) 

COMMON E»P0C1♦POC2»POC3»POC4»POC5,POC6» 

1 PNLl»PNL2»PNL3fPNL4,PNL5»PNL6*RN,C»A. 

2 DENOM 

Z = E(NWB ) - EL 
IF(Z)10.3,3 

3 H = 0.3 

W = 10.0*(Z/3.) 

6 NW = W + 1.0 
Y = 0.0 

DO 206 L = 1»NW 

P(L) = C*(Z - Y)*PIN*EXPF(A*(Y**RN - Z**RN )) 

206 Y = Y + H 
SUM I = 0.0 
NY = W 

DO 207 L = 2»NY 

207 SUM I = SUM r + P(L) 

PNO = H/2.*(P(1) + 2.*SUM I + P(NW)) 

GO TO 2 
10 PNO = 0.0- 
2 CONTINUE 
RETURN 
END 
END 
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APPENDIX lie 


SUBROUTINE GROPE 


TRIPLE INTEGRALS 


SUBROUTINE GROPE ( ELI ♦ EL2tELS»EL4,ELS,EL6♦ 

1 PNl ,PN3AN,P0CN,PNLN,EPU»EPL»EPU1»EPL1» 

2 NLl ,NL2»NL3»NL4»NL5»NL6*NPQ,NPR, 

3 .NP1,NP2»NP3»NP4,NP5»NP6 ) 

DIMENSION E(102 ) fPNl(101)fOENOMI101) 

COMMON E»P0C1♦POC2»POC3»POC4,POC5»POC6, 

1 PNL1,PNL2»PNL3»PNL4,PNL5*PNL6»RN»C»A» 

2 DENOM 

DO 999 J = !♦ 101 
IF(E(J) - EPL)10,3,3 

3 CONTINUE 

NWC =''10.*( (E( J) - EPL)/3.) + 1.0 
EP = E(J) - EPL 
IF(E(J) - EPU)11»4,4 

4 CONTINUE 

NWB = 10.0*(E(J) - EPLl)/3. + 1.0 
IF(E(J) - EPU1)12,5,5 

5 CONTINUE 
PIN = 1.0 

21 EL = ELI - EPLl 
IF(EL)72,14,14 

14 GO TO ( 141,72 ),NL1 

141 CALL PROBE II ( EL,PIN,PN3A1,NWB ) 

GO TO 22 

72 PN3A1 = 0.0 

22 EL = EL2 - EPLl 
IF(EL)73,15,15 

15 GO TO ( 151,73 ),NL2 

151 CALL PROBE II ( EL,PIN,PN3A2,NWB ) 

GO TO 23 

73 PN3A2 = 0.0 

23 EL = EL3 - EPLl 
IF(EL)74,16,16 

16 GO TO ( 161,74 ),NL3 

161 CALL PROBE II ( EL,PIN,PN3A3,NW8 ) 

GO TO 24 

74 PN3A3 = 0.0 

24 EL = EL4 - EPLl 
IF(EL)75,17,17 

17 GO TO (171,75 ),NL4 

171 CALL PROBE II ( EL,PIN,PN3A4,NWB ) 

GO TO 25 

75 PN3A4 = 0.0 

25 EL = ELS - EPLl 
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IF(EL)76»18,18 

18 GO TO ( 181,76 ),NL5 

181 CALL PROBE II ( EL»P I N , PN3A5,NWB ) 

GO TO 26 

76 PN3A5 = 0,0 

26 EL = EL6 - EPLl 
IF(EL) 77,19,19 

19 GO TO (191,77 ),NL6 

191 CALL PROBE II ( EL,PIN ,PN3A6,NWB ) 

GO TO 78 

77 PN3A6 = 0.0 

78 CONTINUE 
GO TO 993 

12 PN3A = 0.0 
GO TO 995 

IF( J - NPQ ) 993,201,202 

201 PN3AP = PN3AN 

202 PN3AN = PN3AP 

993 CONTINUE 

PN3A = PN3AN/( PN3A1 + PN3A2 + PN3A3 + PN3A4 + 
1 PN3A5 + PN3A6‘) 

995 CONTINUE 

31 EL = ELI - EPL 
IF(EL) 172,42,.42 

42 GO TO ( 421,172 ) ,NP1 

421 CALL PROBE II ( EL,PN3A,POC1,NWC ) 

CALL PROBE II (EL,1.0,POBl,NWC . ) 

GO TO 32 

172 POCl = 0.0 
POBl = 0.0 

32 EL = EL2 - EPL 
IF(EL)173,43,43 

43 GO TO ( 431,173 ),NP2 

431 CALL PROBE II ( EL,PN3A,POC2,NWC ) 

CALL PROBE II (EL,1.0,P0B2,NWC ) 

GO TO 33 

173 POC2 = 0.0 
POB2 = 0,0 

33 EL = EL3 - EPL 
IFIEL)174,44,44 

44 GO TO ( 441,174 ),NP3 

441 CALL PROBE II ( EL,PN3A,POC3,NWC ) 

CALL PROBE II (EL,1,0 , POB3,NWC ) 

GO TO 34 

174 POC3 = 0.0 
POB3 = 0.0 

34 EL = EL4 - EPL 
IF(EL)175,45,45 

45 GO TO ( 451,175 ),NP4 

451 CALL PROBE II ( EL,PN3A,POC4,NWC ) 

CALL PROBE II (EL,1.0,POB4,NWC ) 

GO TO 35 

175 P0C4 = 0,0 
P0B4 = 0,0 
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35 EL = EL5 - EPL 
IF(EL)176,46,46 

46 60 TO ( 461»176 ) ,NP5 

461 CALL PROBE II ( EL♦PN3A»P0C5.NWC ) 
CALL PROBE II (EL♦1.0»P0B5,NWC ) 

GO TO 36 

176 P0C5 = 0.0 
P0B5 = 0.0 

36 EL,= EL6 - EPL 
IF(EL)177,47,47 

47 GO TO ( 471,177 ),NP6 

471 CALL PROBE II ( EL,PN3A,POC6,NWC ) 
CALL PROBE II {EL,1.0,P0B6,NWC ) 

GO TO 40 

177 P0C6 = 0.0 
P0B6 = 0.0 

48 CONTINUE 
GO TO 992 

11 PNlAl = 0.0 
GO TO 994 

IF ( J - NPR ) 992,901 ,902 

901 POCP = POCN 

902 POCN = POCP 
992 CONTINUE 

PNlAl = POCN/( POBl + POB2 + POB3 + 
1 P0B4 + POB5 + POB6 ) 

994 CONTINUE 
EL = ELI 


CALL 

PROBE 

I I 

( 

EL,PNlAl,PNL1,J 

) 

EL = 

EL2 





CALL 

PROBE 

I I 

( 

EL,PNlAl,PNL2,J 

) 

EL = 

EL3 



- 


CALL 

PROBE 

I I 

( 

EL,PNlAl,PNL3,J 

) 

EL = 

EL4 





CALL 

PROBE 

I I 

( 

EL,PNlAl,PNL4,J 

) 

EL = 

ELS 





CALL 

PROBE 

I I 

( 

EL,PNlAl,PNL5,J 

) 

EL = 

EL6 





CALL 

PROBE 

I I 

( 

EL,PNlAl,PNL6,J 

) 


PNl(J) = PNLN/OENOM(J) 
GO TO 991 


10 PNKJ) = 0.0 


991 CONTINUE 

POCl 

= 

0.0 

POC2 

= 

0.0 

POC3 

= 

0.0 

P0C4 

= 

0.0 

POC5 

= 

0.0 

POC6 

= 

0.0 

PNLl 


0.0 

PNL2 

= 

0.0 

PNL3 

= 

0.0 

PNL4 

s 

0.0 

PNL5 

X 

0.0 
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PNL6 = 0.0 
PN3A1 = 0.0 
PN3A2 = 0.0 
PN3A3 = 0.0 
PN3A4 = 0.0 
PN3A5 = 0.0 
PN3A6 = 0.0 
POBl = 0.0 
POB2 = 0.0 
POB3 = 0.0 
P0B4 = 0.0 
POB5 = 0.0 
POB6 = 0.0 
999 CONTINUE 
RETURN 
END 
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APPENDIX III 


1. P(E) 

2. P(E) 

3. P(E) 
k- P(E) 

5. P(E) 

6. P(E) 

7. P(E) 

8. P(E) 

9. P(E) 

10. P(E) 

11. P(E) 

12. P(E) 

13. P(E) 
12f. P(E) 


DENSITY FUNCTIONS 
= C (S - E) EXP( A( e’^- S^)) ; 2 = E -Eth 
= C EXP( A E’^) 

= C 

= C S EXP( A ( E^ - Eth^ )) 

= C S EXP( A E^) 

= C S EXP( A ( e’^ - 2^ ) ) 

= C E EXP( A E^ ) 

= C E EXP ( A ( E^ - 2^ ) ) 

= C E EXP ( A ( e’^ - Eth’^ )) 

= C/EXP( A E^) 

= C E EXP( 2(AE)^ ) 

= C E EXF< 2 AE ) 

= C E/EXP( AE^) 

= C(E)^ EXP( AE^ ) 
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APPENDIX IV 


TEST 

In order to test the program, it was postulated that 
if all threshold values were set equal to the same value, 
the output of the program would be constant and the 
results wotild be the input multiplied by that constant. 

In carrying out this test, the density fvinction: 

P(E) = C E EXP( A E^ ) 
was used. All threshold values were set equal to 7.2 
MEV, the lowest threshold value. 

The result is a graph which is the plot of the input 
multiplied by a constant as expected. See Figure 65 . 
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